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Observational Critical QCD 
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A detailed study of correlated scalars, produced in collisions of nuclei and associated with the 
(j-field fluctuations, {8a) 2 =< a 2 >, at the QCD critical point (critical fluctuations), is performed 
on the basis of a critical event generator (Critical Monte-Carlo) developed in our previous work. The 
aim of this analysis is to reveal suitable observables of critical QCD in the multiparticle environment 
of simulated events and select appropriate signatures of the critical point, associated with new and 
strong effects in nuclear collisions. 
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The existence of a critical point in the phase diagram 
of QCD, for nonzero baryonic density, is of fundamental 
significance for our understanding of strong interactions 
and so its experimental verification has become an is- 
sue of high priority Q. For this purpose an extensive 
programme of event-by-event searches for critical fluc- 
tuations in the pion sector has already been started in 
experiments with heavy ions from SPS to RHIC energies 
. In we have emphasized, however, that in order to 
reveal critical density fluctuations in multiparticle envi- 
ronment, one has to reconstruct the cr-sector from pion 
pairs (<7 — > 7r + 7r~) and study correlations of sigmas as a 
function of the invariant mass near the two -pion thresh- 
old. In fact, the QCD critical point, if it exists, communi- 
cates with a zero mass scalar field (cr-field) which at lower 
temperatures (T < T c ) may reach the two-pion threshold 
and decay in very short time scales owing to the fact that 
its coupling to the two-pion system is very strong. Ob- 
viously, the fundamental, underlying pattern of cr-field 
fluctuations, built-up near the critical point by the uni- 
versal critical exponents of QCD 3], is phenomenologi- 
cally within reach if and only if the study of correlated 
sigmas, reconstructed near the two-pion threshold, be- 
comes feasible. In this Letter we perform a detailed fea- 
sibility study of the observables related to the detection 
of the QCD critical point in nuclear collisions. In order to 
proceed we summarize, first, the principles on which the 
behaviour of a critical system of sigmas is based 0, . 

(a) The geometrical structure of the critical system in 
transverse space (after integrating in rapidity) consists of 
cr-clusters with a fractal dimension dp = ^fep leading 

to a power law, < a 2 >~ Ixp^+i, for the cr-field fluctua- 
tions, within each cluster (<5 : isotherm critical exponent) 

(b) In transverse momentum space the c-fluctuations 



obey a power law < a 2 >^ [pj 



2(8-1) 

" i+i 



leading to ob- 



servable intermittent behaviour of factorial moments: 
F 2 (M) ~ (M 2 )ttt where M 2 is the number of 2D cells 
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(c) The density fluctuations of pion pairs (tt + tt~~) with 
invariant mass close to two-pion threshold (2m,r) simu- 
late to a good approximation the sigma-field fluctuations, 
{8a) 2 «< (t 2 >, at the critical point, under the assump- 
tion that the sigma mass reaches the two-pion threshold 
{ra a *~ Im-x) at a freeze-out temperature close to the 
critical value. 

(d) The QCD critical point belongs to the universality 
class of the 3D Ising system in which 8 w 5. 

On the basis of these principles and the fact that criti- 
cal clusters in the above universality class interact weakly 
[5j, one may construct a generator of critical events (Crit- 
ical Monte-Carlo: CMC) containing only sigmas, corre- 
lated according to the above prescription Q . The numer- 
ical experiment is completed by a self-consistent treat- 
ment of cr-reconstruction in the simulated events (after 
the process of decay of sigmas to pions) and a compar- 
ative study between the associated critical correlations- 
fluctuations and the corresponding behaviour of a con- 
ventional Monte-Carlo (HIJING). 

The input parameters of the simulation are the size of 
the system in rapidity (A) and transverse space (R±) as 
well as the proper time scale (r c ) characteristic for the 
formation of the critical system. The CMC generator 
can be used to obtain a large set of critical events. The 
analysis of these events could serve as a guide for the 
search of the critical point in the data of real A + A 
processes in SPS and RHIC experiments. This will be 
shown in a more transparent way in what follows. 

We will investigate two sets of critical events. One 
set corresponding to a rather small system with mean 
charged pion multiplicity per event < n c h >« 40 and an- 
other set describing a larger system (< n c h >~ 120). 
Each data set consists of 10 5 CMC generated events. 
First we will consider the properties of the small system. 
For the values of the input parameters in this case we take 
A = 6 adapted to the SPS energies and R± = 15 fm. 
In order to reveal underlying critical fluctuations, at the 
level of observation, one has first to perform factorial mo- 
ment analysis in small cells of the phase space 0. We 
have chosen transverse momenta for this analysis, the 
reason being that density fluctuations in this plane are 
practically independent of the geometry of the critical 
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FIG. 1: The second factorial moment Fb in transverse mo- 
mentum space of the negative pion sector produced through 
the decay of the critical sigmas in the 10 5 CMC events. For 
comparison we show in the same plot the second moment for 
the sigmas before their decay. 



system Q. Applying factorial moment analysis to the 
transverse momenta (p x ,Py) of the negative pions in the 
sample of the 10 5 critical events we obtain for the second 
moment the behaviour shown in Fig. 1. We observe a 
weak intermittency effect (F2 ~ M 2s2 with s 2 ~ 0.07) 
much smaller than the expected to occur in the critical 
system (s 2 = 2/3). For comparison we show in the same 
plot the corresponding factorial moment for the sigmas, 
before their decay, which follows closely the theoretical 
prediction: s 2 ~ 0.63. The pions analysed in this figure 
originate from the decay of the CMC sigmas using pure 
kinematics and taking into account the isospin properties 
of the decay amplitude. The mass of the decaying sig- 
mas is treated as an almost uniformly distributed random 
variable 0. 

The reason for the suppression of fluctuations in the pi- 
onic sector is the kinematical distortion of the self-similar 
pattern formed in the sigma sector due to the sigma- 
decay. The strength of thi s distortion increases with mo- 
mentum transfer: Q = y/ m 2 — im 2 , and becomes neg- 
ligible near the two-pion threshold (m CT 2m,,). Here 
rn^ = {p 7T + + p n -) where p^± are the four momenta of 
the produced charged pions. 

In this treatment we consider only the visible mode of 
the sigma-decay to charged pions and the reconstruction 
of the sigmas, based on this sector, is simply a search for 
(7r + , 7r~) pairs with invariant mass close to 2m T . The aim 
of this procedure is to establish a self-consistent mecha- 
nism for producing critical sigmas either directly from the 
event generator (CMC) or through reconstruction from 
the final states (7r + , 7r~ ) in the same data set of simulated 
events. The significance of this self-consistent reconstruc- 



tion approach in a search for critical a- fluctuations in real 
events becomes self-evident. In practice this can be done 
by looking, in each recorded event, for (tt + ,tt^) pairs 
fulfilling the criterion A: 

A = {(ir+,ir-)\4ml < {p w + + p^-f < (2m w + e) 2 } (1) 

with *C 1. The momentum of the sigma-particle 
is then obtained as: p a = p v + + p v - . In general the 
pairing will not be one-to-one, allowing for different pos- 
sibilities to pair the charged pions in order to get sig- 
mas. It turns out that the full pairing forming all possi- 
ble pairs (7r + ,7r~) fulfilling Q for a given 7r + , produces 
more accurately the geometric properties of the critical 
sigma sector in the limit e — > 0. Due to finite statis- 
tics this limit can only be taken using data sets consist- 
ing of events with reconstructed sigmas, obtained for a 
sequence of e- values approaching to zero, and then ex- 
trapolate smoothly to e — 0. Thus we have performed 
factorial moment analysis in transverse momentum space 
using values of e in the region 8 — 15 MeV . We have cal- 
culated the second moment F2 for different sets of events 
with reconstructed sigmas, each obtained using a fixed 
e. It exists a wide window in the resolution scale M for 
which F2 ~ (M 2 ) S2 ( e '. Using linear extrapolation the full 
pairing leads to the value S2(0) = 0.70 ± 0.03 very close 
to the theoretically predicted s 2 (0) = §■ In an alterna- 
tive rule for pairing, choosing for each tt + randomly the 
corresponding n~ partner, from all possible ones fulfill- 
ing and repeating the same moment analysis we find 
s 2 (Q) ~ 0.44 ± 0.01 deviating by 35% from the theoret- 
ically expected value. Thus the full pairing reproduces 
succesfully the geometrical properties of the critical sys- 
tem due to the fact that in this way one reconstructs all 
critical sigmas even if it is impossible to identify them. 

The above algorithm besides the reconstruction of the 
real critical sigmas also creates a number of fake sig- 
mas. Some of these are generated by combining pions 
belonging to different critical sigmas and contribute sig- 
nificantly to the strong intermittency effect occuring in 
the reconstructed isoscalar sector. Therefore the inter- 
mittency effect in the factorial moment analysis is a nec- 
essary but not sufficient condition to identify critical, 
strongly correlated isoscalars. There is a need to sub- 
stract the effect of tfie fake sigmas in order to reveal the 
correlated critical isoscalar sector. To achieve this goal 
one has to calculate the pair-correlation density of the 
sigmas defined by: 



C(l,2) =< P2 (l,2)>-<p( m) (l,2)> 



(2) 



where < P2 > is the two-particle density averaged over 
the events. The arguments 1,2 refer to the momentum 
coordinates of two sigma-particles. The second term is 
the two-particle density estimated using mixed events [|| . 
In fact for a given set of events one can construct a set 
of mixed events by shuffling the momentum coordinates 
of the pions in the initial data set. As a consequence the 
(possible) correlations in the original data, including the 



3 



formation of sigmas, are destroyed in the mixed set. In 
general the two-particle density is given by the number 
of particle pairs per event and per phase space cell of 
volume dfl. Therefore eq.|0) can be written as: 

i=l iy ei) i = l 

where i labels the events, N ev is the total number of 
CMC events fullfiling the condition while is 
the total number of mixed events, obtained from the 
shuffling of the tracks in the CMC events, fullfiling QJ. 
The coefficient W is the ratio of numbers of events: 
W = S° . The correlator G(M) is defined as the inte- 
gral of C(l, 2) over a finite phase space cell with volume 
Sil = -Srr, where fl is the total available phase space in 
D-dimensions (here D = 2). Performing the integration 
and normalizing properly we can express the correlator 
in terms of the second order factorial moments 6] . Tak- 
ing into account the fact that the resulting distributions 
of the sigmas depend on the parameter e through 0) we 
finally get: 

( <n {m) >V < ^ 
G(M, e) = F 2 (M, e) - W(e) F 2 (m) (M, e) 

\ Ma ^ J 

. (4) 

where < n a > is the mean multiplicity of sigmas within 
a phase space cell for the CMC events while < ni > is 
the corresponding quantity for the mixed events. These 
two multiplicities are in fact almost equal and therefore 
their ratio appearing in eq. J2J is very close to one for all 
e. The weight factor W has a transparent physical in- 
terpretation: it describes the relative weight of the fake 
sigmas in the reconstructed set. If W(e) w 1 the fake sig- 
mas dominate while for W(e) — > they are suppressed 
relative to the real sigmas. It turns out that the depen- 
dence of W on e is crucial for the determination of the 
critical isoscalar sector. Accordingly, if W ~ the cor- 
relator follows the properties of the second moment Fi 
and the effect of the mixed events is marginal, while for 
W « 1 the correlator becomes negligible and no critical 
effect can be identified. At threshold (e = 0) we ex- 
pect to fully recover the critical isoscalar sector therefore 
W(0) must vanish if real critical sigmas exist in the data 
set of the reconstructed sigmas. This should be the case 
when analyzing the CMC events. In fact, it turns out 
that in order to avoid fake sigmas at threshold, in the re- 
construction procedure, events with sufficient multiplic- 
ity of pion pairs (n a ) must be considered in the analysis 
(n<T > 2). In addition the correlator G(M,e) should be- 
have like ~ M 2 ^ 2 ^ if self-similar patterns of correlated 
sigmas are present in the corresponding phase space. The 
exponent 02 (e) is a new intermittency exponent having 
the property 2 (O) = s 2 (0) as G(M, 0) « F 2 (M,0). Fi- 
nally, the strength of the correlations in the sigma sector 
is given by the magnitude of G{M, e). Therefore in terms 
of the correlator we have three very strict requirements 



in order to reveal the critical isoscalars in a data set of 
charged pions: 

• Dominance of the real sigmas with respect to the 
fake sigmas when approaching the two-pion thresh- 
old (lim e ^ W(e) = 0) 

• Large positive values of the correlator at threshold 
(G(M, 0) 3> 1) at high resolution scales. 

• Correct power-law dependence of the correlator on 
the resolution scale M, at threshold: G(M,0) ~ 
M 1 - 34 

In a conventional hadronic system with no critical 
isoscalar sector at least one of the above properties will 
not be met. 

To test the above ideas we applied the sigma- 
reconstruction algorithm to 10 5 CMC events describing 
a small system as well as to 33000 HIJING events Q 
for the C + C system at the SPS. The mean charged 
pion multiplicity per event in these two systems is practi- 
cally the same. The CMC data set is expected to possess 
a recoverable critical isoscalar sector. On the contrary 
the HIJING simulator describes, in the low p± region, 
conventional hadronic dynamics and therefore no criti- 
cal isoscalars are expected to occur in the corresponding 
data set 0. 

We have first calculated the correlator for these two 
systems choosing n a > 5 for the multiplicity of pion pairs. 
The results are given in Fig. 2 for e-values in the region 
of 8 — 15 M eV, suggested by finite statistics. 

As we can see the CMC-correlator behaves very differ- 
ently from the HIJING-correlator for all values of e. In 
Fig. 3 we also show the weight factor W(e) for both sys- 
tems. It is obvious that in the CMC case W approaches 
zero for e — > while in HIJING W(e) is of order of one 
independently of the e-value. The interpretation of this 
very important result is that in the case of HIJING, i.e. 
in a conventional hadronic system, the reconstructed sig- 
mas are fake for every value of e. 

Finally in Fig. 4 we plot 02(e) for the CMC case. The 
linear extrapolation leads to 02(0) = 0.69 ± 0.04 in very 
good agreement with the theoretically expected value of 
2/3. According to the analysis above we have recovered 
in a self-consistent manner the critical isoscalar sector in 
the CMC data set as a result of a-reconstruction. On the 
contrary the HIJING events do not contain real critical 
sigmas, as expected. 

We also examined the dependence of our results on 
the charged pion multiplicity. We have produced CMC 
events for a larger system with size in rapidity A = 11 
and R± = 25 fm. In this case the total charged pion 
multiplicity is w 120 per event coming from the decay of 
sigmas with mean multiplicity k 96. We found that W(0) 
is very close to zero (linear extrapolation: W(0) — —0.03) 
and 02(0) = 0.664 ± 0.02 in very good agreement with 
the estimations in the smaller critical system (Fig. 4). 
Thus the results concerning the description of the critical 
isoscalar sector are not sensitive to the size of the system, 
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FIG. 2: The correlator G(A1, e) for the reconstructed sigmas in 10 s CMC events and in 33000 HIJING events. We present the 
results for 4 different values of e. 
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FIG. 3: The weight factor W(e) for the CMC and the HIJING 
events. 

as expected. What is changing is the correlation length 
which for a second order phase transition is proportional 
to the size of the system. In fact for the two systems 
considered here the critical correlation length is w 6 fm 
(small system) and « 10 fm (large system) respectively 
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FIG. 4: The slope </>2(e) of the correlator for the CMC system. 


In order to support further our claim that in the 
CMC description the sector of 7r + 7r~ pairs near two-pion 
threshold, has a dinstict dynamical behaviour, associated 
with (T-field fluctuations at the critical point, we have ex- 
amined the behaviour of tt + tt + pairs in the same region 
(e ~ 4 MeV). As expected no effect is found, the pairing 
7r + 7r + has no dynamical origin (W ~ 1) and the corre- 
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FIG. 5: The correlator G(M, e) for (tt + ,tt + ) pairs (stars) as 
well as for sigmas (open triangles) using e = 4 MeV. 



lator G(M, e) of these pairs is strongly suppressed. We 
illustrate this in Fig. 5 where we compare, for e = 4 MeV, 
the correlators of sigmas (7r + 7r~ pairs) and 7r + 7T + pairs 
for the small system (< n c h >~ 40). It must be noted 



that for this value of e the Coulomb correlations are neg- 
ligible (Q 2 = -( Pn+ -v*-? « 50 MeV 0). Similar 
results are found if the system is larger (< n c h >~ 120). 

In conclusion, we have shown that a set of well pre- 
scribed observables (factorial moments, correlators, sta- 
tistical weight of fake sigmas, intermittency exponents) 
associated with the existence of a critical point in quark 
matter, can be established in nuclear collisions. These 
observables belong to the reconstructed cr-sector de- 
scribing massive scalars (7r + 7r~ pairs) near the two-pion 
threshold and their behaviour reveals strong critical ef- 
fects suggested by a-field fluctuations near the criti- 
cal point. We claim that a search for such a critical 
behaviour in heavy ion experiments is feasible within 
the framework of a reconstruction procedure of massive 
scalars, discussed in this work. Since the critical effects 
in this sector are strong and their pattern remains robust 
for systems of different size, our proposal is to study, us- 
ing the above observables, different processes at the SPS 
and RHIC with the aim to scan a substantial area of 
the phase diagram, in a systematic search for the QCD 
critical point in collisions of nuclei. 
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